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Abstract. The onset of the thermonuclear runaway in a Chandrasekhar-mass white dwarf, 
leading to the explosion as a type la supernova, is studied with hydrodynamical simulations. 
We investigate the evolution of temperature fluctuations ("bubbles") in the WD's convective 
core by means of 2D numerical simulations. We show how the occurrence of the thermonu- 
clear runaway depends on various bubble parameters. The relevance of the progenitor's 
composition for the ignition process is also discussed. 
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1. Introduction 

The transition from the hydrostatic to the ex- 
plosive C-burning in a carbon-oxygen white 
dwarf (CO WD), which accretes mass in a bi- 
nary system until it approaches the Chandrase- 
khar mass, marks the start of the explosion of 
this progenitor star as a type la supernova (SN 
la). From a computational point of view, the 
modelling of this phase is extraordinarily chal- 
lenging because of the span of involved length 
and time scales at ignition. 

Nonetheless there are several studies on 
the evolution of the pr ogenitor to the igni- 
tion of SN la explosion. iLesaffre et alJ (l2006h 
study the accretion phase of the progeni- 
tor and find interesting links between bi- 



nary population synthesis and the physi- 
cal conditions of the WD at ignition. On 
shorter time intervals, the last stage before 
the thermonuclear runaway ha s been explored 
(W 

4)£ 

ulatio ns (iHoflich & SteinI I2002L iKuhlen et all 



en explored 
with analytical m odels ( Wooslev et al.l 120041 
IWunsch & Woosley ,2004) and numerical sim- 
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12006*). Though widely different in their numer 
ical approach, these studies agree in identify- 
ing the convective flow prior to the runaway as 
a crucial issue for the ignition process and the 
early stage of the explosion. 

A complementary way to address the 
ignition problem is given by small-scale 
(length scales of the order of 1 km or 
less, whereas the WD diameter is about 
2000 km) investigations of the evolution 
to runaway (e.g. Garcfa-Senz & Bravol I2005L 
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IZingale & Du rsi'2007). Accordin g to the mod- 



Temperature (K) 



els of 



iG arcfa-Se nz & WooslevI (|l995h and 
IWunsch & Woosle^ 20041) . thethermonuclear 
runaway is triggered by temperature perturba- 
tions ("bubbles") on the km-scale, generated in 
the WD core by the turbulent convection. The 
study of the bubble features is therefore a pow- 
erful tool to better understand the ignition pro- 
cess. 

In this contribution, the numerical simula- 
tions (Sect. |2] and [3]l of llapichino et all (l2006l) 
(henceforth 106) are reviewed (we refer to that 
paper for further details). In Sect. |4] we extend 
the discussion of the results to assess the role of 
the WD composition and th e relevance of n ew 
nuclear measurements ( Spillane et al.ll20(j7h . 




2. Numerical setup 

The 2D simulations were performed using the 
FLASH code (Fryxell et al., 2000) . The bub- 
ble is initialised as a temperature perturba- 
tion, hotter and less dense than the back- 
ground, and in pressure equilibrium with the 
surrounding matter The size of the computa- 
tional domain is 5 x 20 km, the geometry is 
Cartesian. The background WD model (pro- 
vided by S. Woosley) is mapped in plane- 
parallel approxim ation, applying the proce- 
dure described by IZingale et alj (|2002|) to en- 
force the hydrostatic equilibrium. A constant 
gravitational acceleration, computed from the 
WD model, is applied to the domain, with the 
gravitational acceleration pointing downwards 
along the y-axis. 

The adaptive mesh refinement (AMR) is 
used, allowing five levels of refinement with an 
effective grid size of [256 x 1024] zones, cor- 
responding to an effective spatial resolution of 
2 X 10^ cm. 

A Helmholtz EOS (Timmes & Swestvl 

12000ft is used. A minimal a-network with 7 iso- 
topes dTimmes et alJl206ol) is implemented to 
follow the hydrostatic C-burning. 

As inferred by the anal ytical models 
dWunsch & Woosleyl 120041: IWoosley et all 
l2004l) . bubbles in the WD core are generated 
with some range of sizes, temperatures and 
distances from the WD centre. A series of 
thirteen calculations was performed over the 




Fig. 1. Temperature plot of the rising bubble at f = 
0.3 s, with the velocity field superimposed, in a run 
with the initial parameters T = 7.7 x 10' K, D = 1 
km,R= 100 km. 

range of expected values for these relevant 
bubble parameters i.e. the initial bubble 
temperature T (7.3 - 7.9 x 10^ K), the initial 
bubble diameter D (0.2 - 5 km) and the initial 
distance from the WD centre R (50 - 150 km). 
A snapshot of the bubble motion is shown 
in Fig. [T| (see 106 for the complete temporal 
sequence). 

3. The physics of buoyant bubbles 

The evolution of the bubble is governed by the 
interplay between nuclear burning and hydro- 
dynamical instabilities. The hydrostatic carbon 
burning determines the nuclear timescale, de- 
fined as the time needed by the temperature 
perturbation to reach the thermonuclear run- 
away (T ~ lO' K), of the order of some sec- 
onds. On the other hand, the bubble is subject 
to the Rayleigh-Taylor instability. It is acceler- 
ated upwards by the eff'ective gravitational ac- 
celeration ^eff('") - g(r) ^p/p, where g is the 
gravitational acceleration at the distance r from 
the WD centre and Ap/p is the density con- 
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Fig. 2. Evolution of the maximum bubble tempera- 
ture in a simulation with the initial parameters as in 



trast between the bubble and the background. 
The vortical motions produced by the bubble's 
buoyant rise, shown in Fig,[Tl lead to its frag- 
mentation and dispersion down to length scales 
in the range 10 - 100 cm, where the heat loss 
by thermal conduction is larger than the nu- 
clear heating. The dispersion timescale is well 
approximated by the buoyant rise time Tbrt - 
Djvb, where D is the bubble diameter and Vb 
is the terminal bubble velocity (of the order of 
1 kms-i). 

The outcome of the bubble evolution is de- 
termined by these competing effects and, quan- 
titatively, by the corresponding timescales. If 
the nuclear timescale is smaller than the buoy- 
ant rise time, the bubble reaches the ther- 
monuclear runaway before being dispersed. 
Conversely, the bubble is dispersed and cools 
down before a thermonuclear flame is ignited. 
The evolution of the bubble temperature in the 
latter case is shown in Fig,|2] where T increases 
until the dispersion prevails. Further numerical 
issues in this process and limitations in our ap- 
proach are discussed in 106, 

4. Results and discussion 

4.1. The parameter study and its 
implications for SN la ignition 

The role of the bubble parameters on the 
physics of temperature perturbations has been 



explored with the performed parameter study. 
The main results are reported here: 

- The bubble diameter D affects the buoyant 
rise time becau se T^n (vh D^^^ 
as well, as in iGoncharovl |2002|) , Larger 
bubbles have longer Tbrt and are therefore 
favoured to go to thermonuclear runaway, 

- The bubble temperature T affects the nu- 
clear timescale, because of the high sensi- 
tivity of the burning to T. Hotter bubbles 
are favoured to go to runaway. However, 
to some extent, also the buoyant rise time 
decreases with increasing temperature, be- 
cause a larger temperature contrast leads to 
an increase of g^s, producing a more effec- 
tive bubble dispersion, 

- The increase of the central distance R in- 
creases ges and thus decreases the buoy- 
ant rise time, because both the bubble tem- 
perature contrast and g(R) increase with 
increasing R. Bubbles at smaller distance 
from the WD centre are favoured to go to 
runaway. 

The evolution of 2D bubbles is a highly 
simplified model for the small-scale ignition, 
nonetheless the results can be used to provide 
some clues on the ignition process in SNe la at 
larger length scales. First, we observe that the 
estimated convective velocity before the run- 
away (50 - 100 km s"') is much larger than 
the bubble buoyant velocity. The features of 
the convective flow in the progenitor prior to 
runaway are therefore needed to model the ig- 
nition properly. Comparing convective veloc- 
ities and timescales for bubble evolution, we 
estimate that the ignition is initiated at a cen- 
tral distance of about 150 km, with a disper- 
sion of the order of 100 km, which agrees well 
with previou s results ( Wooslev et al, I I2004t 
IWunsch & Wooslev 2004) , 

The bubbles are generated in the WD's 
core with some probability distribution func- 
tion of temperature fluctuations. Our bubble 
analysis indicates that the (few) hotter bub- 
bles are more effectively dispersed than the 
(many) milder temperature fluctuations, be- 
cause of their shorter bubble rise time. It sug- 
gests, indirectly, that the multi-point ignition 
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scenario is a viable model for the onset of the 
SN la explosion. 

4.2. Role of Other parameters 

The study o f the diversity in e xplosion simu- 
lations (e.g. iRopke et al.ll2006l) . caused by the 
physical parameters of the WD and linked with 
the observed diversity in SNe la, has received 
comparatively more attention than the explo- 
ration of the diversity in the ignition phase. The 
bubble model described here has intrinsic and 
numerical limitations in its predictive power of 
the ignition properties of SNe la, but it can be 
profitably used as a "probe" of the influence on 
the ignition process of other physical parame- 
ters. 

In iRopke et al] (|2006|) the role of the 
WD composition is explored by varying the 
carbon/oxygen ratio and ^^Ne mass fraction 
(representative of the progenitor's metallicity, 
[Timmes et al. 2003). Dursi & Timmes (2006) 
found that the ignition timescale depends 
strongly on the -^Ne abundance. We run some 
test simulations which confirm their results in 
the framework of the bubble physics. We sug- 
gest that the importance of both the C/O ratio 
and X(^^Ne) for the ignition process should be 
carefully scrutinised in future works. 

Recently Spil lane et al.l ([20071) presented 
new measurements of the cross sections of the 
reactions ^^CC^C, a)^°Ne and i2c('2c,p)23Na 
from £■ = 2.10 to 4.75 MeV. The authors claim 
that the new reaction rates have interesting as- 
trophysical applications. The improvement of 
the experimental determination of these reac- 
tion rates in the considered energy range (or 
at lower energies) could significantly influence 
the hydrostatic carbon burning in the progeni- 
tor and is worth being explored. The resonant 
screeni ng effects should also be further inves- 
tigated dltoh et al.ll2003l) . 
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